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ABSTRACT 

The hrst star clusters and quasars re- 
suhed directly from the growth of linear 
density fluctuations in the early Universe. 
Since they emerged from a well-defined set 
of initial conditions, the interplay between 
observational data from the Next Genera- 
tion Space Telescope (NGST) and theoret- 
ical modeling will advance the current un- 
derstanding of star and quasar formation. 
The first objects had a substantial impact 
on the thermal and chemical state of the 
rest of the Universe as they reionized the 
intergalactic medium and enriched it with 
metals. 

1. INTRODUCTION 

The Universe followed the arrow of time 
in pedagogical order. It started with the 
simplest initial state of almost perfect ho- 
mogeneity and isotropy and ended up with 
systems as complex as intelligent organ- 
isms at present (Fig. 1). So far, stan- 
dard textbooks are unable to fully fol- 
low this pedagogical track because of the 
lack of observational data about the Uni- 
verse in the redshift interval 5 ^ z ^ 10^. 
The COBE satellite probed the microwave 
anisotropics which are remnant from the 
recombination epoch at z ~ 10^ (Ben- 



nett et al. 1996), and existing telescopes 
are limited in their range to redshifts 
z ^ 5. The Next Generation Space Tele- 
scope (NGST) will bridge the above gap in 
our knowledge and image the first sources 
of light that had formed in the Universe. 
With its exceptional nJy sensitivity in the 
l-3.5//m infrared regime, NGST is ide- 
ally suited for probing optical-UV emis- 
sion from sources at redshifts ^ 10, just 
when popular Cold Dark Matter models 
for structure formation predict the first 
baryonic objects to have collapsed (Fig. 
2). 

The first sources are a direct conse- 
quence of the growth of linear density fluc- 
tuations. As such, they emerge from a 
well-defined set of initial conditions and 
the physics of their formation can in prin- 
ciple be coded to the form of a precise com- 
puter simulation. The cosmic initial con- 
ditions for star formation are much sim- 
pler and better defined than those in the 
local interstellar medium; they include a 
prescription for the primordial power spec- 
trum of Gaussian density fluctuations, the 
mean temperature and density of the gas, 
the light element abundances according to 
Big-Bang nucleosynthesis, and the lack of 
dynamically-significant magnetic fields. 
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Figure 1: Milestones in the evolution of the Universe from simplicity to complexity. 
NGST will bridge between the recombination epoch probed by COBE (z ~ 10^) and 
the horizon of current observations {z ^ 5). 
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The initial mass function of the first stars 
is determined by this set of initial con- 
ditions (while subsequent stellar genera- 
tions arc affected by photoionization heat- 
ing and metal enrichment feedbacks). Al- 
though the early evolution of the seed den- 
sity fluctuations can be fully described an- 
alytically, the collapse and fragmentation 
of nonlinear structure must be simulated 
numerically. The first baryonic objects 
connect the simple initial state of the Uni- 
verse to its complex current state, and 
their study offers the key to advancing 
our knowledge on the formation physics of 
stars and massive black holes. 

2. THE REIONIZATION EPOCH 

The collapse redshifts of baryonic ob- 
jects in a standard CDM cosmology is 
shown in Figure 2. As seen from this 
figure, more than a tenth of all baryons 
could have assembled into virialized ob- 
jects by a redhift z ~ 10. Since nuclear 
fusion yields ~ 7 MeV per baryon and 
black hole accretion might provide even 
more energy, and since hydrogen ioniza- 
tion requires 13.6 eV, it follows that even 
the collapse of only a negligible fraction 
( ^ 10^^) of all baryons into stars or black 
holes could ionize the rest of the Universe. 
A detailed calculation of the expected star 
and quasar populations in the popular A- 
CDM cosmology, implies that the Universe 
is typically reionized at z ~ 10 (Haiman & 
Loeb 1998a). 

Reionization is defined as the time when 
the volume filling factor of ionized hy- 
drogen (HII) approached a value close to 
unity. The transition to this state was 
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Figure 2: Collapse redshift, Zcoii, for cold 
dark matter (dashed lines) and baryons 
(solid lines) in spheres of various baryonic 
masses, Mb, and initial overdensities. The 
overdensities are in units of the rms ampli- 
tude of fluctuations (t(M) for a standard 
CDM power-spectrum with a^h-\ = 0.67. 
The collapse of the baryons is delayed rel- 
ative to the dark matter due to gas pres- 
sure. The curves were obtained by fol- 
lowing the motion of the baryonic and 
dark matter shells with a spherically sym- 
metric, Lagrangian hydrodynamics code 
(Haiman & Loeb 1997). 
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DETERMINING THE REIONIZATION REDSHIFT 
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Figure 3: Sketch of the expected spectrum of a source at a redshift Zg sHghtly above 
the reionization redshift Zreion- The transmitted fluxes due to HII bubbles in the pre- 
reionization era and the Lya forest in the post-reionization era are exaggerated for 
illustration. 
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almost sudden, because the ionizing 
sources were typically separated by a dis- 
tance which is much smaller than the Hub- 
ble length so that the HII regions they 
had produced overlaped on a timescale 
much shorter than the Hubble time. Is 
it possible to identify the reionization red- 
shift, Zrcioii) from the spectrum of high red- 
shift sources? The most distinct feature 
in the spectrum of a source at Zs > -Zreion 
would be the Gunn-Peterson absorption 
trough due to the neutral intergalactic 
medium that fills the Universe prior to 
reionization. Figure 3 provides a sketch 
of the spectrum of a source with (1 -|- 
Zs) > 1.18 X (1 + 

2^rcion) which Contains 
some transmitted flux between the Gunn- 
Peterson troughs due to Lja and LyP ab- 
sorption. This transmitted flux is sup- 
pressed by the residual Lya forest in the 
post-reionization era. The possibility of 
identifying Zreion from the damping wing 
of the Gunn-Peterson trough (Miralda- 
Escude 1997) suffers from potential con- 
fusion with damped Lya absorption along 
the line of sight and from ambiguities due 
to peculiar velocities and the proximity ef- 
fect. An alternative method makes use of 
the superposition of the spectra of many 
sources (Haiman & Loeb 1998b). In the 
absence of the Lya forest this superpo- 
sition should result in the sawtooth tem- 
plate spectrum shown in Figure 4 [with the 
horizontal axis redshifted by (1 -l-^reion) to 
the observer's frame]. The cumulative ra- 
diation produced by the first sources might 
also get reprocessed through intergalactic 
dust and appear as a diffuse infrared back- 
ground (Haiman &; Loeb 1998a). 



The first stars have also enriched the 
intergalactic medium with metals. The 
observed disparity (Lu et al. 1998) be- 
tween the metallicity of the Lya forest 
( ^ 1Q~^Zq) systems and the damped 
Lya absorbers (~ Q.\Zq) indicates that 
some of the metals are retained within 
the potential wells in which they are pro- 
duced. Nevertheless, the nonzero mctalic- 
ity found in the intergalactic gas implies 
that winds driven by supernovae, massive 
stars, gamma-ray bursts (Loeb & Perna 
1998), or merger events (Gnedin &; Os- 
triker 1997) expelled some of the enriched 
gas out of galaxies. A typical Type H 
supernova releases ^ 1 eV in thermal 
energy per ~ 1Q~^Zq in metals. The 
virial temperature of an object of total 
mass Mtot which formed at a redshift Zf 
is ~ 1 eV (Mtot/lO^Me)'/^ [(1 + Zf)/10]. 
Hence, gas could 

be expelled out a potential well only if its 
metallicity exceeds the minimum value of 
- 10-3^0 (Mtot/IO^M©)'/-' [(1 + Zf)/1Q]. 
Moreover, in order for mixing to occur ac- 
cross the Hubble velocity difference be- 
tween star clusters, one needs to sup- 
ply the gas with a kinetic energy of ~ 
1 eV X [(1 + z)/l{)f (d/10 kpc)^ where d 
is the average distance between neighbor- 
ing star clusters. Since objects of mass 
M collect their mass from a radius r ~ 
10 kpc (Mtot/10^Me)i/''^[(l + z)/10]-i and 
since at high redshifts d ^ r, considerable 
mixing will occur only if the metallicity in- 
side the objects is increased well above the 
level of 10^'^ Zq. The early stars are ex- 
pected to behave as coUisionless particles 
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Figure 4: The combined spectrum of many pre-reionization sources, ignoring absorp- 
tion by the Lya forest after reionization. The upper panel shows that absorption by 
neutral hydrogen and helium suppresses the flux above 13.6cV up to the keV range. 
The lower panel shows the sawtooth-modulation due to line absorption below 13.6eV 
(from Haiman, Rees, & Loeb 1997). 
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and eventually populate the halos of 
present-day massive galaxies. Therefore, 
direct observations of the abundance pat- 
terns in the halo stellar population of the 
Milky- Way could uncover the early enrich- 
ment history of the first star clusters and 
the intergalactic medium. 

3. DETECTION WITH NGST 
How many sources will NGST see? Fig- 
ure 5 shows the predicted number of 
quasars and star clusters expected per field 
of view of NGST (Haiman & Loeb 1998b). 
In this calculation, the star formation ef- 
ficiency was calibrated based on the in- 
ferred metallicity range of the Lya for- 
est (Songaila & Cowie 1996; Tytler et 
al. 1995) while the characteristic quasar 
lightcurve was calibrated in Eddington 
units so as to fit the observed luminos- 
ity function of bright quasars at z 2- 
4. Both populations of sources were ex- 
trapolated to high redshifts and low lu- 
minosities using the Press-Schechter for- 
malism (for more details, see Haiman &: 
Loeb 1998a). Typically, there are of or- 
der tens of sources at redshifts z > W per 
field of view of NGST. The lack of point 
source detection in the Hubble Deep Field 
is consistent with a low-mass cutoff for lu- 
minous matter in halos with circular ve- 
locities ^ 50-75 km s~^, due to photoion- 
ization heating (Haiman, Madau, & Loeb 
1998) . The redshift of early sources can be 
easily identified photometrically based on 
their Lya trough. Figure 5 demonstrates 
that NGST will play a dominant role in 
exploring the reionization epoch and in 
bridging between the initial and current 
states of the Universe. 
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Figure 5: Predicted number counts per 
5' X 5' field of view per logarithmic flux in- 
terval in the NGST wavelength range of 1- 
3.5/im. The numbers of quasars and star 
clusters were calculated for a ACDM cos- 
mology with {QM,^A,^h,h,agh-i,n) = 
(0.35, 0.65, 0.04, 0.65, 0.87, 0.96). The low- 
est mass scale of virialized baryonic ob- 
jects was chosen consistently with the pho- 
toionization feedback due to the UV back- 
ground. The star formation efficiency was 
calibrated so as to bracket the possible val- 
ues for the average metallicity of the Uni- 
verse at z ~ 3, namely between 10~^Zq 
and W~'^Zq. The thick lines, labeled 
"10", correspond to objects located at red- 
shifts z > 10, and the thin lines, labeled 
"5", correspond to objects with z > 5. 
The upper labels on the horizontal axis 
correspond to Johnson I magnitude (from 
Haiman & Loeb 1998b). 
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